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ABSTRACT 

A stellar population synthesis model, suitable for comparison with Giant Extragalac- 
tic Hn Regions (GEHRs), is constructed incorporating the recent developments in mod- 
elling stellar evolution by Maeder and co-workers and stellar atmospheres by Kurucz. A 
number of quantities suitable for comparison with broad band data of GEHRs in visible 
and near infrared parts of the spectrum are synthesized in addition to the hydrogen 
and helium ionizing photon production rates at solar metallicities, for three scenarios 
of star formation — (i) Instantaneous burst (IB) (ii) Continuous star formation (CSF) 
and (iii) Two bursts of star formation, with the older burst rich in red supergiants. For 
IB case, evolution of colors shows three distinct phases — an initial steady blue phase, 
followed by a red bump (5-15 Myr) and another steady phase with colors intermediate 
to the earlier two phases. CSF colors asymptotically reach peak values at ~ 10 Myr, 
never reaching the reddest IB colors. Ionizing photon production rate falls off by an 
order of magnitude in 6 Myr for IB, where as it almost remains constant for CSF model. 
Two-burst models with burst separations ~ 10 Myr have properties of both IB and 
CSF, simultaneously producing the red IB colors and high ionizing photon rate, making 
such regions easily distinguishable using optical observations. Flat IMFs result in bluest 
colors when the massive stars are on the main sequence and reddest colors during the 
red supergiant phase of the evolving massive stars. Errors on the computed quantities 
due to the statistical uncertainties inherent in the process of star formation become neg- 
ligible for cluster masses in excess of 10 5 M & . Our GEHR spectra in the range 200 nm 
to 3 /mi are found to be in good agreement with the computations of Mas-Hesse and 
Kunth (1991). 
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1. Introduction 

External galaxies offer excellent opportunities for the 
study of star formation under a variety of physical 
conditions. Many of the external galaxies experience 
a much higher level of star formation activity com- 
pared to our own Galaxy. Regions of current star for- 
mation are concentrated around the nucleus in star- 
burst and Sersic-Pastoriza galaxies (Balzano 1983, 
Sersic & Pastoriza 1965), while they are spread over 
the entire disks in blue compact dwarfs, also known 
as Hn galaxies (Terlevich et al. 1991). Giant ex- 
tragalactic H II regions scattered over the disks of ir- 
regular galaxies and arms of spiral galaxies are the 
other regions associated with intense on-going star 
formation activity. At the presently achievable spa- 
tial resolutions, individual stars in most of the active 
star-forming regions cannot be resolved or counted, 
owing to their greater distances. Thus, information 
regarding the stellar population and the star forma- 
tion history in these regions has to be obtained from 
the integrated properties of all the embedded stars. 
Such a study would involve a comparison of synthe- 
sized properties of star-forming regions with observ- 
ables. Most of the early efforts in this direction con- 
centrated on synthesizing integrated colors of galaxies 
(Larson & Tinsley 1978, Searle et al. 1973), which are 
dominated by old low mass stars (age ~ 10 10 yr). 
Bruzual & Chariot (1993), Buzzoni (1989) synthe- 
size integrated spectra of galaxies, which again sample 
populations averaged over lifetimes of galaxies. The 
effect of active star formation lasts at most for a few 
tens of millions of years, with significant variations 
over period of fractions of million year (Myr). Thus 
it is important to have time steps of the order of frac- 
tions of Myr to obtain a complete picture of the early 
evolution of an active star-forming region. Extensive 
grid of models covering a variety of observables in se- 
lected bands from ultraviolet (UV) to near infrared 
(NIR) have been synthesized by Rieke et al. (1980), 
Lequeux et al. (1981), Melnick et al. (1985), Olof- 
sson (1989), Mas-Hesse & Kunth (1991) and more 
recently by Leitherer & Hcckman (1995), all suitable 
for young star forming complexes. 

In all the above models a giant star-forming com- 
plex is assumed to consist of a star cluster with an as- 
sociated nebulosity around it. Ionizing photons from 
hot OB stars of the star cluster are the source of ion- 
ization of the nebula. Majority of the emission lines 
from the complex are produced in the gaseous neb- 
ula, while the UV-optical continuum has contribu- 
tion from both stars and ionized gas. The dust par- 
ticles embedded inside the H n regions absorb a frac- 
tion of ionizing photons, re-emitting them in far in- 



frared wavelengths. The models from different groups 
differ mainly in the input data, and in the slightly 
different synthesis techniques used. There are two 
main approaches adopted in the literature for syn- 
thesizing observable quantities from star-forming re- 
gions. Bruzual & Chariot (1993) and Mas-Hesse & 
Kunth (1991) among others have used observed stellar 
quantities as input data base, while Buzzoni (1989) 
and Melnick et al. (1985) have used theoretical atmo- 
spheric models of Kurucz (1979) as input. Although 
the former method is devoid of uncertainties in the 
model atmospheres, it suffers from the incompleteness 
in the stellar input data base, especially for the fast 
evolving evolutionary stages. In the present work, we 
follow the second method, which gives the freedom to 
study the active star-forming regions over a variety 
of physical conditions. The main aim of the present 
model is to synthesize the spectra and selected observ- 
able quantities in the optical and NIR regions, incor- 
porating the latest developments in modelling stellar 
evolution and stellar atmospheres. Specifically the 
computations are done to compare the observables 
in the optical bands obtained through CCD imaging. 
The most important of these quantities is the ratio of 
H/3 to blue band luminosity. This ratio resembles the 
H/3 equivalent width, but unlike equivalent widths, it 
can be obtained accurately through imaging observa- 
tions. It should be noted that H/3 equivalent widths 
are normally obtained through slit spectroscopy, with 
slits often small compared to the spatial extent of the 
giant star-forming complexes. This necessitates the 
computations of quantities suitable for imaging ob- 
servations. 

We are considering the phase of the starburst in 
which the most recent one or two bursts dominate 
the emission in all the bands from UV to NIR. This 
phase corresponds to ages of the order of a few times 
10 7 yrs. We consider three scenarios of star formation, 
and evolve the cluster up to 20 Myr after the epoch 
of star formation. In the first scenario, all the clus- 
ter stars are assumed to have formed simultaneously 
(Instantaneous Burst). Stars are assumed to be form- 
ing continuously (Continuous Star Formation) in the 
second scenario. Third scenario we considered has a 
second burst occurring in the star formation complex 
before the death of low to intermediate mass stars 
from the earlier burst (multiple bursts). 

In section 2, we define the model star- forming com- 
plex, and describe the methods used for synthesizing 
the stellar quantities. The section also includes the 
comparison with observed stellar quantities. The re- 
sults of the computations for the three scenarios dis- 
cussed above is presented in detail in section 3, while 
the effect of statistical uncertainties inherent in the 



2 



process of star formation are discussed in section 4. 
The computed quantities are compared with the re- 
sults from earlier work in section 5. Main results are 
summarized in section 6. 

In paper I (Mayya 1994) of this series we reported 
photometric data in BVRHa bands for around 180 
Hn complexes in nine galaxies. In the third and con- 
cluding paper in this series we will derive parameters 
related to star formation by comparing the observa- 
tional data of Paper I with the synthesis results of 
this paper. 

2. Description of Synthesis Technique 

The present model assumes the star forming com- 
plex to be a star cluster surrounded by an ionized 
medium. Its spectrum is thus obtained by adding 
the spectrum of an ionized nebula to a pure cluster 
spectrum. This section deals with the description of 
the synthesis technique used in computing these two 
spectra. 

2.1. Definition of a Cluster 

A star cluster in the model is defined by the fol- 
lowing set of parameters : 

1. Initial Mass Function (IMF) of stars defined by 

(a) Slope (a) 

(b) Upper cut-off (ra„) 

(c) Lower cut-off (mi) 

2. Age (t) 

3. Metallicity (Z) 

4. Total Mass (M T ) 

The interesting range of above parameters are, m u = 30- 
120M Q , a =1.0-3.5, t =0-20 Myr, Z =0.001-0.04 
(Z = 0.02 corresponds to solar metallicity, Z & ). 
Lower cut-off mass mi can have any value greater than 
O.1M . Mi is one of the most difficult parameters to 
be evaluated, because of the dominance of flux from 
higher mass stars in any wavelength region of obser- 
vation. However, a significant fraction of the total 
mass of star-forming regions is locked up in low mass 
stars and hence the total mass as well as star forma- 
tion rates are critically dependent on mj . Based on a 
number of observable quantities, Rieke et al. (1980) 
found mi to be in the range 2-3 Mq in starburst re- 
gions. In the years that followed, there has been a 
number of studies favoring mi values in the range 1- 
10M©(see Scalo 1987 for a discussion). We assume 
a value of 1 M & for mi throughout the present work. 
Effect of changing mi on the computed quantities is 
discussed in sec. 3.1. Following formulations are used 
in the computations of the spectral energy distribu- 



tion and luminosity in any band from a star cluster. 

m u 

L x (imf;t;Z) = A J l x (m;t; Z)m~ a dm (1) 

mi 

and 

L B (imf;t;Z)= A J l B (m;t; Z)m~ Q dm (2) 

mi 

with 

M T = A J m 1 - a dm (3) 

mi 

l x (m;t;Z) = irf x (m;t;Z)4irR 2 (m;t;Z) (4) 

where f\(m; t; Z) is the intensity emitted from a star 
of mass m at an evolutionary stage t at wavelength A. 
R(m; t; Z) is the radius of the star, A, the normalizing 
constant. Z is the abundance of the gas from which 
the stars are formed. The stellar luminosity in any 
given band B with a normalized response curve i?s(A) 
is computed as, 

l B (m;t;Z)= J l\(m;t; Z)R B (\)d\. (5) 

band 

Lyman continuum photon rate is obtained by inte- 
grating the spectrum between 22.8 nm and 91.2 nm. 
Thus, 

91.2 

n L (m;t-Z)= J l x (m;t; Z)^-d\, (6) 

22.8 

with h and c being Planck's constant and speed of 
light respectively. The above equations assume the 
cluster to be sufficiently massive so that masses can 
be considered to be continuously populated, instead 
of discretization into mass bins or spectral types. Un- 
der these conditions, numerical integration can be 
performed rather than summing over all the mass 
points. In sec. 4, we evaluate the minimum mass of 
the cluster above which numerical integration rather 
than summation over mass is justifiable. 

2.2. Input Stellar Data 

The computations require a library of stellar spec- 
tra for a wide range of masses at various evolution- 
ary stages at the metallicities of interest. We used 
the theoretically computed stellar spectra as input 
dataset to our models. These spectra are synthe- 
sized by combining the stellar evolutionary models 
from Maeder and co-workers and stellar atmospheric 
models from Kurucz. We have used the evolutionary 
tracks from Schaller et al. (1992), which is their most 
recent model for Z—Zq. These grids differ from the 
earlier grids from the same group (Maeder & Meynet 
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1989, Maeder 1990) mainly in having incorporated the 
new radiative opacities from Rogers & Iglesias (1992). 
The models assume convective core overshooting dis- 
tance of 0.20 times the pressure scale height, and a 
metallicity dependent mass loss rate (M ~ Z - 5 ). For 
masses greater than 20M Q , the authors have tabu- 
lated evolutionary results for two values of the mass- 
loss rate, one being twice the other. The latter mod- 
els are computed with the intention of reproducing 
the properties of the low luminosity Wolf-Rayet stars. 
As we are not aiming to count the number of Wolf- 
Rayet stars as the cluster evolves, the two mass-loss 
schemes are expected to give identical results. We 
have used the set with lower mass-loss rates, which 
is more standard of the two. Effective temperatures 
for Wolf-Rayet stars in these models have been com- 
puted taking into account the extended atmospheres 
as described in Kudritzki et al. (1989). Massive star 
evolution is carried out until the end of core carbon 
burning in these models. 

The atmospheric models of Kurucz (1979) also 
have been recently updated (Kurucz 1992). Apart 
from including many more layers than before to im- 
prove the numerical accuracy of the computations, 
the new models have a resolution of 1-2 nm in UV 
and optical regions of the spectrum. New computa- 
tions are available at 1221 wavelength points from 9 
nm to 160 /im. The effective temperature T c ff range 
of 3500-50000 has been covered for logarithm surface 
gravity(logg) values of 0.0 to 5.0. Kurucz grids are 
adequate for all stages during the evolution of stars 
in the mass range 1-120, except in the following situ- 
ations. Effective temperature of main sequence stars 
more massive than 80M Q (earlier than 04) exceed 
50000 K. These very hot models do not differ much 
from blackbody and hence we used blackbody equa- 
tions when ever the temperature exceeded 50000 K. 
At the low temperature end, during the fast evolving 
phase of red supcrgiants, the temperature drops below 
3500 K, for initial masses m < 20 at ages > 10 Myr. 
During this phase the log g values also cross the lower 
limit of in Kurucz grids. In these cases we have used 
the 3500 K Kurucz models with \ogg — 0. This ap- 
proximation will not affect our computations for ages 
< 10 Myr. Beyond this age, small effects might be 
seen in the NIR. 

2.3. Computational details 

Using Schallcr et al. (1992) evolutionary models, 
we have computed \ogg and T e g at the 51 tabulated 
evolutionary stages for Z = Z Q and 1 < m < 120. 
The spectrum for a given log g and T c g is computed 
using Kurucz (1992) atmospheric models at Z — Zq, 
interpolating the spectrum linearly in T e g and log g 



when necessary. The absolute luminosities l\ are com- 
puted from these spectra by using the mass-radius re- 
lationship as derived from the tables of Schaller et al. 
(see eq. 4) . This results in a spectrum sampled at all 
the 1221 Kurucz wavelength points. 

It can be seen from equations 1 and 2 that com- 
putations of cluster quantities require a knowledge of 
l\ (or ^b) at every mass point chosen for integration. 
Schaller et al. tabulate the evolutionary grids for 21 
mass points in the mass range 0.8-120M Q , which is 
inadequate for a complete integration over the chosen 
mass limits. Thus it is desirable that the luminosities 
1a are computed at sufficient number of mass points 
to represent every phase during the evolution of the 
cluster. This necessitated the generation of tracks for 
stellar masses intermediate to the tabulated masses. 
After considerable experimentation we found it nec- 
essary to extend the 21 mass Schallcr et al. grids 
to 151 mass grids. Of these 102 mass points lie be- 
tween 10-40-Mq, which is the most important mass 
range for evolutions upto 20 Myr. Sampling interval 
is smoothly increased from 0.2 to O.5M in this mass 
range. For smaller number of grids the evolutionary 
results are found to be sensitive to the exact choice of 
the mass points. The interpolation scheme adopted in 
generating the evolutionary tracks for any intermedi- 
ate mass is as follows. For a given mass, tracks for two 
adjacent masses from Schaller et al. are identified. A 
track is characterized by the mass, effective temper- 
ature (T c ff) and bolometric luminosity (Zboi) and is 
sampled at 51 evolutionary stages. New tracks for in- 
termediate masses are generated by linear interpola- 
tions in log(T c ff)-mass and log(^boi)-mass planes at all 
the 51 evolutionary stages. The tracks so obtained are 
sampled at 0.1 Myr interval up to 20 Myr by linear in- 
terpolations in log(T c ff)-age and log(^boi)-age planes. 
The chosen time step is suitable enough to represent 
all the short-term stellar evolutionary phases. The 
band luminosities /^(m; t; Z) arc further interpolated 
at mass intervals of 0.05 dex for m < 40M Q and 0.2 
dex for higher masses to carry out the numerical in- 
tegrations. Interpolations are done linearly in log(^s) 
- log(m) plane. 

The interpolation scheme described above is adopted 
after visually verifying the linearity of the interpo- 
lated quantities in the plane of interpolation. 

2.4. Synthesized Quantities and Comparison 
with Observations 

Spectral energy distribution covering wavelengths 
in UV, optical and NIR regions for any given stel- 
lar mass is one of the basic quantities synthesized in 
our model. Luminosities in BVRIJHKLM bands 



4 



and the number of hydrogen and helium ionizing pho- 
tons are synthesized from these spectra using the cqs 
in sec. 2.1. The response curves needed are taken 
from Bessell (1990) for BVRI bands and Bessell & 
Brett (1988) for JHKLM bands. Band luminosities 
are transformed into colors using the zeropoints de- 
rived from the Vega spectrum of Kurucz (1992), as- 
suming Vega to have colors 0.00 mag in all bands. 
Zeropoints are also derived using the observed Vega 
spectrum for BVRI bands and model spectrum of 
Dreiling & Bell (1980) for JHKLM bands. The ze- 
ropoints derived from these two methods agree within 
0.01 mag. Lyman continuum photon rate is obtained 
by integrating model spectra shortward of 91.2 nm as 
given in Eq. 6. Similarly the helium ionizing luminos- 
ity is computed by integrating spectra shortward of 
50.4 nm. The absolute luminosity scale in the model 
is tied to Schaller et al. (1992) models. This approach 
does not require the use of the highly uncertain T e g- 
Bolometric Correction(BC) relation, to set the lumi- 
nosity scale of the model. On the other hand, we 
derive the BC over the region of interest in T e ff and 
log g in order to compare with the observed values. 

Before discussing the computations of cluster quan- 
tities, we compare the computed model colors and 
BCs with observed values for different spectral types 
and luminosity classes. The most homogenous set of 
stellar photometric data in BVRIJHKL bands are 
due to Johnson (1966). Flower (1977) has updated 
the T e ff.B — V:BC scales making use of the most re- 
cent photometry available at that time, and found 
Johnson (1966) scale to be in general agreement in 
most cases, with B — V colors redder by around 0.12 
mag for giants. For comparison with the model results 
we used B — V and BC relations with T e ff from Flower 
and the rest of the photometry from Johnson. John- 
son V — R colors are transformed into Cousins system 
using the transformation equations given by Bessell 
(1979) (after flipping the mis-printed sign of the ad- 
ditive constant for the redder color range: Bessell, 
private communication) . There is more data available 
in the literature for the evolving cool and hot phases 
of the stars (Humphreys & McElroy 1984; Elias et 
al. 1985), which are in general agreement with the 
compilations of Johnson (1966) and Flower (1977). 

Lyman continuum photon production rate per unit 
surface area (ni yman ) from our computations is com- 
pared with that from Mas-Hesse & Kunth (1991) in 
Fig. 1(a). Other computed quantities such as bolo- 
metric correction (BC), B—V, V — R and V— K colors 
are compared with the observed values in Fig. 1(b)— 
(e) as a function of T c ff . In each of these figures, the 
observed values are denoted by filled circles for main 
sequence stars and inverted triangles for supergiants. 



Main sequence model values are shown by the solid 
lines, whereas the dashed line represents the super- 
giant values as derived from our computations. It can 
be seen that on the whole, there is good agreement be- 
tween our model computations and observations. Our 
Teffii? — V scale for supergiants is around 0.03 mag 
bluer than that of Flower (1977), with good agree- 
ment for main sequence stars. The differences are of 
the order of 0.1 mag for red supergiants below 10000 
K, where the color-T e ff relation is very steep. 5-10% 
errors in T c ff for a given color can easily explain these 
differences. V — R and V — K show similar trends, 
with the model colors redder by as much as 0.2 and 
0.5 mag respectively for cool red supergiants. The ob- 
served BC for main sequence and supergiants lie on 
either side of our computations, with the differences 
of the order of 0.2 mag for T e ff greater than 10000 
K. For the coolest red supergiants the model BCs are 
higher (less negative) by as much as 0.5 mag. As in 
the case of B — V, 5-10% shift of model T e gf scale 
towards the cooler side for a given color can explain 
the observed differences in V — R 1 V — K and BCs. 
Such shifts in T e ff are not unreasonable considering 
the different ways in which T e g can be defined for 
extended atmospheres of stars, especially supergiants 
(Baschck et al. 1991). The Lyman continuum rates 
agree within 0.05 dex for hot main sequence stars, 
which mainly contribute to the total ionizing photons 
in a cluster. The observed differences at lower temper- 
atures where the Lyman continuum steeply decreases 
with T e g are thus not much significant. 

In our model, the observational quantities are pa- 
rameterized in terms of mass and hence the absolute 
luminosities depend on the mass-luminosity relation. 
As mentioned earlier, our scale is based on stellar evo- 
lutionary models of Schaller et al. (1992), which are 
found to be in good agreement with the values tabu- 
lated by Lang (1992) for main sequence stars. Similar 
comparison for supergiants could not be performed 
because of the uncertainty in the mass estimates for 
the model supergiants. However model derived T c ff- 
luminosity relation agrees with the observed positions 
of red supergiants in the Hertzsprung-Russel diagram 
(Humphreys & McElroy 1984). Total Lyman contin- 
uum photon production rate for high mass main se- 
quence stars agrees within 0.1 dex with the computa- 
tions of Mas-Hesse & Kunth (1991), Leithcrcr (1990) 
and Panagia (1973). Our values begin to get system- 
atically different (lower) at masses lower than 4OM , 
with differences reaching a factor of 2 for 2OM . The 
ionizing photon rate is lower by about a factor of ten 
for these lower masses compared to those from the 
hottest stars. 

The basic stellar quantities and computed lumi- 
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nosities for the selected stellar masses at zero age are 
given in Table 1(a) for solar metallicity. Table 1(b) 
lists the same quantities during the red and blue su- 
pergiant phases of the stars. The latter table contains 
quantities which are obtained by taking the weighted 
averages over the supergiant phases. The quantities 
tabulated in the two tables are: 

Mass, M ms : Zero age main sequence (ZAMS) mass 
of the star 

M sg : Mass of the star during its supergiant phase 
T : Time spent in main sequence phase (Tab. la) and 
supergiant phases (Tab. lb) 

T e ff, BC, rad, Mboi : Effective temperature, Bolo- 
metric Correction, radius and bolometric magnitude 
of stars 

: Ratio of stellar mass to the luminosity in the 
visual band, both expressed in solar units. M ms is 
used as mass of the supergiants rather than the highly 
mass- loss dependent M sg . 

Nl : Number of ionizing photons emitted from the 
star per second 

/n e : Fraction of lyman continuum photons (A < 
91.2 nm) available for the ionization of helium also 
(A < 50.4 nm) 

: Ratio of the expected H/3 luminosity from the 
nebula surrounding the star to the blue band lumi- 
nosity. See sec. 3.1 for an exact definition of this 
quantity, -f- , /h c and Nl are meaningful only for 
very hot stare (T off > 25000K) 
colors : Optical and NIR broad band colors. 



2.5. Computations of Nebular Flux 

Ultraviolet photons from hot stars in the cluster 
ionize the surrounding gas leaving hydrogen, helium, 
oxygen, nitrogen, carbon and other astrophysically 
abundant elements in different ionized states. For ion- 
ization bounded nebulae, all of these ionizing photons 
are used up by hydrogen, and at ionization equilib- 
rium, every ionization is followed by recombinations 
to one of the levels in the hydrogen atom (Osterbrock 
1989). In thermal equilibrium, the energy input to 
the nebula in the form of kinetic energy of the elec- 
trons is lost due to collisions and radiative recombi- 
nations, leaving the nebula at a fixed temperature. 
This heated gas along with the ions in various ex- 
cited states gives rise to a spectrum characteristic of 
Hn regions. The basic physics governing the emis- 
sion mechanism from these regions is well understood 
and hence its theoretical spectrum can be synthesized 
from only a few input parameters such as temper- 
ature, density, abundance of elements and absolute 
flux in one of the Balmer lines. In the next two sub- 
sections, we describe the method we have followed in 



synthesizing the continuum and line spectrum of the 
nebula. 

2.5.1. Estimation of the Gaseous Continuum 

The main emission mechanisms contributing to 
continuum emission from a gas are free-free, free- 
bound, and two-photon (2q) emissions from hydro- 
gen and helium. The net emission coefficient j c g can 
be expressed as a sum of these emission coefficients, 
weighted by the relative abundance of the particular 
ionized state. 

7eff - 7(HI)+7(2q)+^^7(HeI)+ n(HC ^ +) 7 (HeII), 



n(H) 



n(H) 



(7) 

where n(H) is the total number density of hydro- 
gen atoms in both ionized and neutral states. In 
our computations we assume a helium abundance of 
10% by number all of which is in the singly ionized 
state. Individual 7's as computed and tabulated by 
Brown & Mathews (1970) are used here for A < lfim. 
For longer wavelengths computation of Ferland (1980) 
have been used. The continuum of the gas spectrum 
in erg s _1 A -1 is related to j e g by, 



Lneb 
A 



7eff 




1 



E H 0<^(HO,T e )^ 



-■H/3- 



(8) 



Er/3 is the energy of the H/3 photon, a^ /3 (H° ,T e ) is 
the H/3 recombination coefficient and Lh/3 is the H/3 
line luminosity in erg s _1 . T e is the temperature of 
the nebula at thermal equilibrium. 

2.5.2. Estimation of the Gaseous Emission line Strengths 

Osterbrock (1974, 1989) has compiled the emission 
line strengths of various transitions of hydrogen and 
helium relative to H/3 line strength for typical densi- 
ties and temperatures of Hn regions. The strength of 
emission lines from other elements however are depen- 
dent on the relative abundance of the elements, the 
state of ionization, the temperature structure of the 
nebula etc. and hence their computations require de- 
tailed modelling. Such models have been constructed 
by Dopita & Evans (1986), Rubin (1985), McCall et 
al. (1985) among others. From observations of a large 
sample of extragalactic H 11 regions McCall et al. have 
shown that the majority of the emission line ratios 
can be parameterized by a single parameter, namely 

I[OII]+I[OIII] 



I(H/3) 



The basic physics behind this can be 
summarized as follows. The cooling of the nebula 
is mainly controlled by the radiative losses in oxy- 
gen ions, which constitute more than 45% (Lcqueux 
et al. 1979) of metallic mass. This gives rise to an 
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anti-correlation between oxygen abundance and neb- 
ular temperature. Most of the oxygen in GEHRs is in 
singly ([Oil]) and doubly ([Oiii]) ionized states, both 
of which have bright emission lines in the optical re- 
gion ([On]A3726 and [Om] A5007). Hence, oxygen 
abundance can be quite accurately estimated using 
one or both of these lines. The observed correlations 
of line ratios of other elements with oxygen abun- 
dance suggests that the relative abundance of other 
elements such as C, N and S do not vary much from re- 
gion to region. Photoionization models are successful 
in reproducing the observed correlation between line 
ratios and oxygen abundance, enabling the prediction 
of line ratios from a knowledge of oxygen abundance 
alone. 

Our interest here is to estimate the contribution 
of emission line fluxes to BVR band fluxes. We con- 
sidered all the bright lines of N, O and S apart from 
hydrogen lines as computed by McCall et al. (1985). 
Their tables 15 and 12 help us to estimate the oxy- 
gen abundance from the strengths of [O n] A3726 and 
[O in] A5007 lines, which in turn are used to estimate 
the strength of various lines from other elements. 

The nebular spectrum is then obtained by adding 
selected emission lines to the continuum spectrum. 
Nebular spectrum when added to the cluster spec- 
trum results in a spectrum expected from a star- 
forming complex. 

3. Star Formation Scenarios 

We consider three possible scenarios of star for- 
mation which might be at work in GEHRs. In the 
simplest case referred as Instantaneous Burst (IB), 
all the stars are assumed to have formed over a short 
time scale (coeval) compared to the lifetime of even 
the most massive stars. The other extreme case is to 
assume the star formation to be continuous over a pe- 
riod of time with a star formation rate (SFR), which 
is either a constant (CSFR) or exponentially decreas- 
ing (ESFR) function of time. Thirdly we consider 
an intermediate scenario between the two mentioned. 
We model the effect of two instantaneous bursts sep- 
arated over time scales comparable to the lifetime of 
the stars produced in the burst. Thus at any given 
time, there will be stars from two generations. We 
refer to this model as Multiple Bursts (MB). The re- 
sults from these three models are described separately 
in the following sections. 

3.1. Instantaneous Burst of Star Formation 

The IB of star formation is fully described by the 
equations in sec. 2.1. All the quantities synthesized 
for individual stars are also synthesized for clusters. 



The dependence of the computed quantities on the 
IMF parameters for a zero-age cluster is presented in 
Tab. 2(a). Model numbers are given in the first col- 
umn. The IMF parameters, namely the upper cut-off 
mass and slope are given in columns 2 and 3. Ra- 
tios of the total mass to bolometric and visual 
(j^-) luminosities are given in columns 4 and 5 in so- 
lar units. Lyman continuum photon production rate 
(Nl) per unit mass of the cluster is given in column 6. 
Fraction of the hydrogen ionizing photons which are 
also capable of ionizing helium are given in column 7. 
Lyman continuum photon rate Nl is also expressed 
in units of expected H/3 luminosity <f) from a nebula 
at 10000 K and defined as, 

<j) = 4.78 x 10~ 13 iV L . (9) 

(j) in the above equation has dimensions of ergs -1 M©~ 
and Nl in phs _1 M© _1 . Column 8 contains , 
which is the ratio of the expected H/3 to blue band 
luminosity. H/3 equivalent width is directly related 
to as can be seen from Fig. 2. The relationship 
is derived by evolving a cluster up to 20 Myr, with 
the dots on the curves separated by 0.5 Myr interval. 
Column 8 contains this ratio. Columns 9-14 contain 
optical and NIR colors of the cluster (nebular contri- 
bution not taken into account). H/3 and Ha equiva- 
lent widths in angstrom are given in column 15 and 
16 respectively, assuming a gas temperature of 10000 
K. 

The effect of changing lower cut-off masses on the 
computed quantities is illustrated in Tab. 2(b), ex- 
pressed as differences from an IMF with mi — 1M© 
in logarithm or magnitude units. The tabulated val- 
ues correspond to an upper cut-off mass of 60M©, but 
have identical values for 30 < m u < 120M©. Tabu- 
lated quantities in columns 4-16 are identical to those 
in Tab. 1, except in column 7, which represents the 
factor by which the mass of the cluster decreases due 
to the change in mi from a value of 1 M©. It can be 
seen that the quantities which are most affected by 
changing the lower cut-off mass are log(^r), log(^) 
and Nl- Changing the lower cut-off has only the ef- 
fect of changing the total mass of the cluster with 
negligible change in the luminosities, as can be in- 
ferred by comparing columns 4-6 with column 7. An 
IMF with a lower cut-off of 1M© contains 4-6 times 
more mass compared to an IMF with m; = 10M© for 
IMF slope = 2.5 and 32-38 times more for slope = 3.5 
for m u in the range 120 to 30M Q . Among colors 
V — K is most sensitive to changes in m ; due to non- 
negligible contribution to K band luminosity from low 
mass stars. However contribution from stars less mas- 
sive than 1M© is negligible to any band luminosity 
and hence colors will remain the same. Quantities in 
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columns 4-6 can be reproduced for a different mi by 
adding the entries in Tab. 2(a) to the logarithm of the 
quantity, 

/m =^£^- do) 

A real giant star-forming complex differs from its def- 
inition in the model in the following ways, 
f . Dust can be present inside star- forming complexes, 
which absorb and hence reduce the total number of 
photons available for ionization. Throughout their 
computations Mas-Hesse & Kunth (1991) have as- 
sumed that 30% of the ionizing photons are absorbed 
by dust, based on the work of Bclfort et al. (1987). 
However Mezger (1978) has shown this fraction to be 
highly uncertain having values anywhere between 0.2 
to 0.7 depending on the local conditions. Because 
of the large uncertainty in the value, we prefer to 
present all the model values without correcting for 
dust absorption. We leave it as a free parameter to 
be assumed or determined while comparing with ob- 
servations. 

2. Reddening due to foreground interstellar dust and 
nebular continuum emission in different bands make 
the observed colors different from the values com- 
puted for a pure cluster. Both these effects can be 
handled in a better way using available information 
on the physical conditions of the individual regions 
while comparing with observations. Thus we tabu- 
late the pure cluster values instead of taking into the 
effects of interstellar extinction and nebular emission 
in an average way. 

3. Coevally formed stars take different amounts of 
time to reach the ZAMS. For example a star of 100M Q 
takes only 1.9 x 10 4 yr, where as a star of 1M takes 
75 Myr (Ezer & Cameron 1967, Beech & Mitalas 
1994). Schaller et al. (1992) models we have used 
considers the stellar evolution from ZAMS stage on- 
wards, thus our early stage cluster evolution might 
be affected. Using the pre-main sequence (PMS) evo- 
lutionary tracks of Ezer & Cameron (1967), we esti- 
mated the PMS contribution relative to the ZAMS lu- 
minosities in various bands. PMS luminosity of even 
the low mass stars remains within a factor of 2-3 of 
the ZAMS value for more than 80% of the PMS evo- 
lutionary stage. On the other hand, in the earlier 
stages, cluster luminosity is dominated by the lumi- 
nous massive main-sequence stars, whose bolometric 
luminosity is 10 5 10 6 times that of 1 M star, ren- 
dering the inclusion of PMS stages unnecessary. Thus 
computations of cluster quantities based on Schaller 
et al. models are valid even at stages where low mass 
stars are yet to reach the ZAMS. 

More massive a star is, earlier it exhausts the nu- 
clear fuel and hence the evolution of the IB model of 



SF is characterized by the selective removal of massive 
luminous stars from the cluster. As the massive stars 
wander around in the H-R diagram during their post 
main-sequence evolution, they affect the integrated 
properties of the cluster, even when their number is 
small. We have evolved the cluster at steps of 0.1 Myr 
interval upto 20 Myr, when the highest surviving mass 
is ~ 11M Q . A detailed discussion of the these results 
will be presented in sec. 3.4 along with the other two 
scenarios of star formation considered. 

3.2. Continuous Star Formation 

In the previous section we considered the evolution 
of a cluster containing coeval stars formed in an in- 
stantaneous burst. Now consider the case where stars 
are formed continuously resulting in stars of differ- 
ent ages in the same cluster. This scenario, known as 
continuous star formation (CSF), is characterized by 
a star formation rate (SFR), the rate at which new 
stars are added to the cluster. We have considered 
two kinds of SFRs, a constant rate of star formation 
(CSFR) and an exponentially decreasing star forma- 
tion rate (ESFR). 

Computationally CSF scenario is realized by sum- 
ming luminosities over all surviving masses, which is 
mathematically represented as, 

t m u 

A j j l B (m;t';Z)m~ a R(t-t') dm dt' 

Lf(imf;t;Z)= ; (11) 

J R(t')dt' 
o 

with 

R(t') = R for CSFR (12) 

and 

R(t') = R exp(-£) for ESFR (13) 

R(t') represents the rate of star formation at an epoch 
t' with Rq and t 0} the SFR at t = and the e- folding 
time of SFR respectively. The denominator in eq. 11 
is a measure of total mass in all the stars and A is the 
normalization constant defined in eq. 3. 

3.3. Multiple Bursts of Star Formation 

We now consider the third scenario of star for- 
mation, in which a second burst of star formation 
takes place before the death of intermediate-low mass 
stars from the earlier burst. In several galactic OB 
associations, neighbouring regions of star formation 
are found to differ in age by ~10 Myr, which is 
thought to be evidence for sequential star formation 
(Elmegreen 1992). The energetic events such as stel- 
lar winds, supernova activity, typical of massive stars, 
might be responsible for the sequential star formation 
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(Elmegreen & Lada 1977). In distant star forming 
regions, the resolution may not be good enough to 
resolve the populations spatially, and hence one has 
to infer the presence of more than one burst based on 
the integrated light. In this study we aim at provid- 
ing observational plots where regions with two bursts 
of star formation can be identified. 

Computations are done by adding the fluxes from 
two populations of equal strength, one younger than 
6 Myr and the other older than 6 Myr. Each popula- 
tions is assumed to be independently evolving like an 
IB with identical IMFs. 

3.4. Discussion of Evolutionary Results 

As discussed in section 2.1, an IB cluster is defined 
by its age, metallicity and total mass apart from the 
three IMF parameters. CSF and MB scenarios have 
additional parameters defining them. Thus the total 
parameter space involved is large for full presentation 
of the results. It should also be noted that Mas-Hesse 
& Kunth (1991) and more recently Leitherer & Heck- 
man (1995) have presented results for considerable 
section of the parameter space. In this study, we con- 
centrate on comparison of the three scenarios of star 
formation for typical values of the parameters, instead 
of presenting results over the entire parameter space. 
We discuss the IB scenario in greater detail for two 
reasons — (i) The evolutionary trends in this case can 
be directly associated with the evolution of stars and 
(ii) The results of IB scenario can be used to explain 
the composite models such as CSF and MB. 

As has been pointed out earlier, we have fixed the 
metallicities at Z =Zq and normalized the quanti- 
ties to cluster mass Mr = 1M Q . An IMF with 
mi = l,m u = 60, and slope = 2.5, which has slope 
close to Salpeter's value of 2.35 for the solar neigh- 
bourhood, is chosen as a representative IMF in our 
computations. We refer this IMF as "standard" in 
rest of the article. The effects of different IMF pa- 
rameters arc inferred by computing quantities for two 
extreme IMFs — (i) mi = l,m u = 30, slope = 3.5, 
and (ii) mi = l,m u = 120, slope = 1.0. The first 
one is heavily biased against high mass stars while 
the second one is rich in massive stars. Comparison 
of the results for the three IMFs is done in sec. 3.6. 

We depict the contribution from different masses 
to the cluster spectrum as a function of age of the 
IB in Figs 3(a)-(f). The mass range 1-6OM , for 
the standard IMF is divided into four sub-intervals, 
namely (i) low mass stars (1 < m < 5), (ii) inter- 
mediate mass stars (5 < m < 15), (iii) high mass 
stars (15 < m < 40) and (iv) very high mass stars 
(40 < m < 60). The quantities chosen are (a) Bolo- 



metric luminosity, (b) Lyman continuum photon rate, 
(c)-(f) B, V, R and K band luminosities respectively. 
The ordinates are expressed in units of erg s _1 M Q _1 
in all figures except (b), where it is in phs -1 M _1 . 
The dotted, dashed and dash-dotted lines represent 
the contribution from stars in 40-60M Q , 15-40M Q 
and 5~15Mq range respectively. The contribution 
from the 1~5M bin always lies below the plotted 
range. The solid line in the plots is the sum of contri- 
bution from all the stars, which is plotted separately 
for four of the optical and NIR bands in Fig. 3(g). 
The bands are marked on the curves. Also marked in 
Fig. 3(g) are the beginning and end of the supergiant 
phases for the selected masses. The following conclu- 
sions can be drawn from these figures. 

1. Early cluster evolution is dominated by the evo- 
lution of stars more massive than 15M Q , with lower 
mass stars beginning to contribute substantially at 
around 10 Myr. 

2. All quantities except Nl show an initial increase 
in luminosity peaking at an age between 5-8 Myr be- 
fore the decline starts. The rate of increase is mainly 
controlled by the evolving massive stars (> 40 Mq). 
Peak luminosities are reached earlier at shorter wave- 
lengths compared to that at longer wavelengths, as 
the IB cluster evolves. Lyman continuum photon rate 
remains steady for the first 3 Myr, decreasing mono- 
tonically then onwards. 

3. The initial increase in luminosity in optical and 
NIR bands is due to the combined effects of the in- 
crease in bolometric luminosity and the decrease in 
the effective temperature of massive stars as they 
evolve off the main-sequence. The increase in lumi- 
nosity is larger at longer wavelengths (e.g. K band). 

4. Bolometric luminosity and ionizing photon rates 
are heavily dependent on the massive luminous stars 
and hence decrease rapidly once these stars die at 
around 5 Myr. On the other hand the luminosities in 
optical and NIR bands decrease only marginally after 
the death of massive stars. 

5. The peak luminosities in optical and NIR bands 
during the evolution of the cluster are due to the red 
supergiant phase of the 15-40 M Q stars. The red su- 
pergiant contribution to total luminosity as compared 
to that of the low mass main sequence stars is increas- 
ingly significant at longer wavelengths. For example 
99% of the K band luminosity is contributed by the 
red supergiants at around 10 Myr. 

6. The only regime where masses less than 5M Q show 
their presence is in the evolution of optical luminosi- 
ties beyond 12 Myr. Even here their contribution is 
less than 20%. 

7. Small wiggles in the diagram especially beyond 
15 Myr are not genuine, but are due to the limit a- 
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tions of the computations. 

The evolutionary results for the three scenarios of 
star formation are presented in Tabs 3, 4 and 5. Refer 
Tab. 2(a) for the explanation of the column headers in 
these tables. Tabulated quantities are smoothed over 
a period of 0.5 Myr to suppress the short term fluc- 
tuations on the computed quantities. For continuous 
star formation model, constant SFR is assumed in the 
computations of the tabulated quantities. Columns 1 
and 2 of Tab. 5 give the ages of the young and old 
bursts for the two-burst models. The data are pre- 
sented for the younger bursts at ages 0, 3.5, 5.0 and 
6.0 Myr and older bursts of equal strength but evolved 
to ages between 6 and 14 Myr. Computational re- 
sults for only the standard IMF are presented in the 
printed form. These Tabs 3-5, containing the results 
for the three IMFs discussed above arc presented in 
AAS CD-ROM series. 

The Evolution of the most important observable 
quantities is shown graphically in Figs 4(a)— (d) for 
the standard IMF. The solid lines depict the quan- 
tities smoothed over 0.5 Myr with the dotted lines 
representing the un-smoothed quantities at intervals 
of 0.1 Myr for IB scenario. Results for CSF scenario 
are plotted as dashed and dot-dash lines correspond- 
ing to exponentially decreasing SFR with to = 10 Myr 
and constant SFR respectively. Quantities plotted are 
log(^ ), B-V, V-R and V-K in Figs 4(a), (b), 

(c) and (d) respectively. ) drops by two orders 

of magnitude in 10 Myr in the case of IB scenario. In 
contrast this quantity falls by only a factor two during 
the same time period for CSFR and by a factor of e 
for ESFR. For IB scenario, follows the evolution 
of Nl, which is controlled by the massive stars. The 
continuous creation of massive stars in CSF models 
keeps the almost steady. 

Three distinct phases can be identified in the evo- 
lution of optical and NIR colors for IB. Colors remain 
constant at their zero age values for the first 3 Myr, 
during which period main sequence stars dominate 
the cluster evolution. The second phase is character- 
ized by a red bump during 5-15 Myr. This phase 
corresponds to the appearance of red supergiants in 
the cluster. Colors again become steady at around 
16 Myr, at a value mid-way between the blue colors 
of main sequence phase and the average colors of red 
supergiant phase. During this phase surviving stars 
have masses <12 Mq and hence the cluster lacks lu- 
minous red supergiants. The color evolution is deter- 
mined by the slowly evolving intermediate mass stars 
and hence colors are steadier. Larger amplitude of 
the dotted lines during this phase represents the limi- 



tations of the computations. Increasing the thickness 
of the mass grid for stars less massive than 12M is 
required to reduce these fluctuations. 

Evolution of the colors in the red supergiant phase 
requires a detailed discussion. B — V attains the red- 
dest value of 0.38 at 7 Myr, drops down to a value 
of 0.12 at 9 Myr, shows two small peaks there after 
before entering the third phase. V — R shows the 
first peak at 8 Myr reaching peak color of 0.4, which 
drops down to 0.29 at 9 Myr. Notice that the drop is 
stronger for B — V compared to V — R and V — K . 
Infact V — K color becomes reddest at 11.5 Myr. The 
position of the color peaks depends on the detailed 
evolution of massive red supergiants, which, when 
present dominate the luminosity in optical and NIR 
bands. Fig. 3(g) is useful in understanding the color 
behaviour. 

Colors for CSF remain steady at zero-age values 
upto 3-4 Myr, increasing to values of 0.1, 0.2 and 2.0 
for B — V, V — R and V — K colors respectively. Colors 
remain steady at these values after about 10 Myr. 
CSF colors remain bluer than the IB colors during the 
luminous red supergiant phase of cluster evolution. 
Difference between CSFR and ESFR models is very 
small compared to their differences from IB models. 

3.5. Cluster Evolution in Observational Plane 

So far we studied the behaviour of the observable 
quantities as a function of age of the cluster. Now 
we discuss the cluster evolution in log(-^ ) vs color 
(B — V, V — R and V — K) plane. These are shown in 
Figs 5(a)-(c) for the standard IMF. Observed quan- 
tities can be directly plotted on these diagrams, and 
hence the plots are very useful for diagnosis of star 
formation. Dotted line denotes the locus of IB model 
evolved upto 13 Myr, with the thicker dots placed 
every 0.5 Myr interval. CSF model with a constant 
SFR is shown by the dashed line, covering 20 Myr of 
evolution. The solid lines correspond to models with 
two bursts of star formation. Four sets of models are 
plotted. Age of the younger population is fixed for 
each set at 0, 3.5, 5 and 6 Myr from top to bottom 
respectively. These ages are chosen to represent most 
prominent evolutionary phases of the cluster and the 
corresponding positions are marked by asterix on the 
IB model. The age of the older population varies from 
6 to 14 Myr along a given set, with thick dots denot- 
ing 0.5 Myr intervals. 

CSF models occupy only a narrow strip in the ob- 
servational plane even after evolution upto 20 Myr. 
On the other hand 2-burst models cover a much larger 
range in log(-^ ) and colors. It is interesting to note 
that these composite models simultaneously produce 



10 



high values and red colors, when the age of the 
younger population is < 3 Myr. Younger population 
controls the ordinate, while the colors are controlled 
by the red supergiants of the older population. This 
phase of cluster evolution can be easily distinguished 
from IB and CSF models. 

3.6. Effect of IMF on cluster evolution 

All the discussions in the previous sections were 
centered on the intermediate IMF. In this section we 
compare those results with two extreme IMFs in order 
to judge the effect of IMF differences on observational 
quantities. The high mass deficient IMF has parame- 
ters mi = 1, m u = 30 and a = 3.5, where as the high 
mass rich IMF has parameters mj = 1, m u = 120 and 
a = 1. 

Evolution of log(^ ), B — V, V — R and V — K 
for IB with the three IMFs discussed above are shown 
in Figs 6(a)- (d). Solid line uses the standard IMF, 
where as the dotted and dashed lines correspond to 
high mass enriched and high mass deficient IMFs re- 
spectively. The initial steady phase in and colors 
prolongs up to 5.5 Myr for the high mass deficient 
IMF. is lower by an order of magnitude during 
this steady phase compared to the standard IMF due 
to the reduction in the number of stars contribut- 
ing to ionizing photons. The red bump in the color 
evolution is less pronounced for this high mass trun- 
cated IMF. Note that the colors are redder during the 
steady phase and bluer in the red supergiant phase as 
compared to the standard IMF values. 

Evolution with high mass enriched IMF is inter- 
esting in several respects. An examination of Fig. 3, 
which is for the standard IMF, is useful in understand- 
ing the details of the evolution. The initial steady 
phase in is almost absent where as it is shortened 
for colors. This is due to the domination of opti- 
cal and NIR luminosities by stars more massive than 
40 Mq. The extent of domination is more for the high 
mass enriched IMF compared to standard IMF. Dur- 
ing the early evolution optical and NIR luminosities 
increase rapidly where as Nl remains steady result- 
ing in the decrease of ■ Due to the shorter lifetime 

° Lb 

of high mass stars, the cluster gets blue and red su- 
pergiants earlier on during its evolution, leading to 
shortening of the initial steady phase in colors. All 
the stars more massive than 60 Mq cease to exist after 
3.5 Myr and hence beyond this age, the high mass en- 
riched IMF differs from standard IMF only in having a 
flatter slope. Evolution of is almost independent 
of IMFs after 5 Myr. The shape of color evolution 
during the red supergiant phase (5-15 Myr) is only 
weakly dependent on IMF parameters, but the abso- 



lute value of colors during this phase is strongly IMF 
dependent. B — V, V — R and V — K colors are redder 
by 0.25, 0.15 and 0.4 mag respectively in comparison 
to the standard IMF. Redder colors are due to an in- 
crease of high luminosity red supergiants with respect 
to low luminosity red supergiants and main sequence 
stars, as a result of the flatter slope. In summary, 
increasing the fraction of high mass stars makes the 
cluster colors bluer in the early phase (< 3 Myr) and 
redder during the red supergiant phase (5-15 Myr). 

Effect of IMFs for CSF models are shown in Figs 7(a)- 
(d). Line types and notations are as given for Fig. 6. 
(-^- ) has a smooth evolution with high mass en- 
riched IMF always having higher value of . Mas- 
sive star enriched IMFs have bluer colors during the 
first 5 Myr and redder colors beyond that compared 
to other IMFs. Reddest colors reached by the CSF 
models are not as red as that of IB models (0.1 vs 
0.65, 0.2 vs 0.55 and 2.0 vs 3.6 for B - V, V - R and 
V — K respectively). 

Two-burst model results for the three IMFs are 
plotted in the observable planes in Figs 8(a)- (c). The 
line types used to distinguish IMFs are as for the pre- 
vious figures. There are four sets of models for each 
IMF. Burst parameters are similar to that for Fig. 5. 
Age of the younger population is marked on each 
track. The high mass enriched IMF spans a consider- 
able range in the observable plane. The three IMFs 
are clearly distinguishable for younger burst ages less 
than 3.5 Myr. Beyond this age, they are separated 
only in the optical colors. 

4. Statistical Effects on the Derived Quanti- 
ties 

Following the normal approach, we have assumed 
the stellar masses to be continuous and have numer- 
ically integrated the Eqs 1 and 2 over mass. Two 
clusters with identical IMFs can have different distri- 
bution of the stellar masses, leading to different band 
luminosities, summed over all the stars. The mean 
of the luminosities over several clusters with identi- 
cal IMFs, is however expected to be the same as that 
obtained using the integration approach. The depar- 
ture from the mean value is dependent on the cluster 
mass Mr, approaching zero for massive clusters. We 
investigate here the minimum mass required for an as- 
sumed IMF, such that statistical fluctuations in star 
formation does not affect the observable quantities. 

We aim to estimate the errors on the computed 
quantities due to the statistical selection of stellar 
masses in a real star forming complex, as a function 
of the total mass of the star cluster. This requires a 
knowledge of the masses of all the stars in the cluster. 
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Monte Carlo method is used for simulating real star 
clusters satisfying the IMF with mi = 1, m u = 120 
and a = 2.5. The IMF is normalized such that at 
least one star earlier than 04 (90 < m < 120) ex- 
ists in the cluster. The mass of the cluster containing 
only one such star is the minimum mass of the cluster 
for the assumed IMF and hence is referred as M m i n 
. The corresponding number of stars in the cluster 
is N min . For the chosen IMF, M min = 6641M Q and 
N m - m = 2434. Cluster masses (Mr) and total num- 
ber of stars (Nt) are defined in integer units of M m i n 
and N m i n respectively. Thus Mr = nxM m - m and 
Nt — nxN m i n . Physically n denotes the number of 
stars earlier than 04. 

For a given mass Mr, Nt stars with masses be- 
tween 1-120M Q are generated weighted by the IMF 
slope. The slope of the IMF for the generated mass 
string is rc-dctcrmincd and confirmed to be within 
0.01 of the input value of 2.5. Luminosity in a given 
band B from such a mass string is computed by 
adding the band luminosities from all the stars. 

/ 1 \ Nt 

L B (imf;t;Z;M T )= f — J £ Z B (m i; t; Z) (14) 

where rrn is the mass of the i th star generated dur- 
ing the simulation. Lb is computed at Z —Zq for 
t =0-20 Myr at intervals of 0.1 Myr. The process of 
generating Nt stars and computing Lbs are repeated 
for 100 trials, thus simulating 100 clusters with the 
same IMF and total cluster mass. The error due to 
the statistical selection of stars is estimated as rms 
noise of the 100 trials. 

The above exercise is carried out for n = 3, 10 and 
50. The mean values of the computed quantities over 
100 trials agree well with the values obtained from 
the integration method. Rms errors on the represen- 
tative luminosities and colors are plotted in Figs. 9 
and 10 for n = 3, 10 and 50. Errors on the bolomet- 
ric luminosity, Lyman continuum photon production 
rate, expressed in terms of H/3 flux, B and K band 
luminosities are plotted in Fig. 9(a)-(d) respectively. 
Errors on , B — V, V — R and V — K colors 

Lb 1 7 

are plotted in Fig. 10(a)-(d). Decrease of rms errors 
with increasing cluster mass (or increasing n) can be 
easily inferred from these plots. The errors on col- 
ors are larger during the red supergiant phase of the 
cluster, as compared to other phases. For n = 50 
(M T = 3.3 x 10 5 M Q ) the errors on log(^ ) are less 
than 0.03 dex throughout the 20 Myr evolution. Er- 
rors remain below 0.03 mag for B — V and V — R col- 
ors during the same period. V — K errors marginally 
exceed 0.1 mag for ages greater than 16 Myr. The 
errors on log(-^- ) and colors are larger by a fac- 



tor of 2-3 for clusters 5 times less massive (n = 10). 
However these errors are within the various photo- 
metric error sources involved during the photometry 
of GEHRs (Mayya 1994). Hence for embedded clus- 
ters more massive than ~ 10 5 M Q or more luminous 
than 10 39 erg s _1 in Ha, the errors due the statisti- 
cal selection of masses are negligible. For comparison, 
30 Doradus nebula in the Large Magellanic Cloud has 
Ha luminosity of 1.5 x 10 40 erg s _1 (Mathis et al. 
1985). 

5. Comparison with Other Models 

Currently a number of models exist, which syn- 
thesize the cluster parameters at various wavelength 
bands. Among these, the most suitable model for 
comparison with our results is from Mas-Hesse & 
Kunth (1991). We compare our monochromatic spec- 
trum of a cluster-l- nebula with their table 6(a). Com- 
parisons are made at all the tabulated ages (2, 5, 7.5, 
10 and 20 Myr) for IMFs with to, = 2, m u = 120 
and a = 2.0. In their model, Mas-Hesse & Kunth 
(1991) have assumed that only 70% of the available 
Lyman continuum photons are used up for ioniza- 
tion. We have taken this into account in our com- 
parison. The results of the comparisons are given in 
Figs ll(a)-(e). Spectra integrated over the entire star 
forming complex are shown in these Figs, between 
70 nm and 3.0 /im at solar abundance, for IB model at 
ages 2, 5, 7.5, 10 and 20 Myr. Solid line shows cluster 
spectrum from our computations, whereas the dotted 
lines show the total cluster spectrum as computed by 
Mas-Hesse & Kunth (1991). The rms of the differ- 
ences between the two spectra over the wavelength 
range 136.2 nm to 2.2 fim arc 0.08, 0.10, 0.18, 0.20 
and 0.08 dex for 2, 5, 7.5, 10 and 20 Myr respectively. 
Hence the agreement is generally good over the en- 
tire wavelength range at all comparison ages except 
at 7.5 and 10 Myr, at which time cluster is in its 
red supergiant phase. Considering that the spectral 
shape is rapidly changing during the red supergiant 
phase and that Mas-Hesse & Kunth derived spectrum 
from available broad band photometry of red super- 
giants, the observed differences during this phase are 
not un-reasonable. 

6. Summary 

We have constructed an evolutionary population 
synthesis model incorporating the recently available 
stellar evolutionary models and model stellar atmo- 
spheres. We synthesize a number of observable quan- 
tities in the optical and near infrared part of the spec- 
trum, suitable for comparison with the observed prop- 
erties of giant star forming complexes. The ratio of 
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H/3 to B band luminosity, B - V, V - R and V - K 
colors have been studied as a function of cluster evolu- 
tion for Instantaneous Burst and continuous star for- 
mation models. Computations performed for a two- 
burst model are also presented with the younger burst 
rich in massive stars capable of ionization, and the 
older burst rich in red supergiants. Color evolution 
shows an initial steady phase for both instantaneous 
burst and continuous star formation models, followed 
by a red bump between 5-15 Myr for instantaneous 
burst. In the observable plane, two-burst models oc- 
cupy quite distinct regions compared to instantaneous 
burst and continuous star formation models. Depen- 
dence of computed quantities on the IMF parameters 
are also studied. Flat IMF is found to produce the 
bluest colors when the cluster is young, and the red- 
dest colors when red supergiants start appearing in 
the cluster. The effect of the statistical fluctuations 
inherent to the process of star formation on the com- 
puted quantities are studied using Monte Carlo simu- 
lations and letting the cluster evolve upto 20 Myr. For 
GEHRs with Ha luminosities exceeding 10 39 erg s^ 1 , 
the statistical effects are found to be negligible. The 
spectra from our models are compared with the com- 
putations of Mas-Hesse & Kunth (1991) at 2, 5, 7.5, 
10 and 20 Myr, and the agreement is found to be 
reasonably good. 
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Fig. 2. — Equivalent width of H/3 emission line is 
plotted against the ratio of H/3 to blue band lumi- 
nosity ) in log-log plot. Age is the parameter 
which is varying along the locus, with successive dots 
separated by 0.5 Myr. The relationship is not very 
sensitive to IMF parameters or the scenario of star 
formation and hence is useful in converting one form 
of data to another. 



Fig. 3. — Contribution to total luminosities from 
stars in four mass bins are plotted as a function of 
cluster age for an instantaneous burst with IMF pa- 
rameters mi — l,m u — 60, and slope = 2.5. The 
chosen mass bins are: 40 < m < 60 (dotted line), 
15 < m < 40 (dashed line), 5 < m < 15 (dot-dashed), 
m expressed in units of M@. The contribution from 
the fourth bin (1 < m < 5) lies below the plotted 
range. The quantities plotted are (a) bolometric lu- 
minosity, (b) lyman continuum photon rate, (c) B 
(d) V (e) R and (f) K band luminosities. The solid 
lines correspond to the summed contribution from all 
the stars which are again plotted in Fig. (g) for a 
direct comparison of band luminosities. The begin- 
ning and the end of the supergiant phases for selected 
masses are marked on this plot. Ordinates are in units 
of erg s _1 Mq _1 for all plots except (b) where it is 
ph s~ 1 M Q ~ 1 . See section 3.4 for further details. 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v3.0. 



Fig. 1. — Lyman continuum photon rate at the 
stellar surface (ni yman ), bolometric correction (BC), 
B — V,V — R, and V — K colors of stars as computed 
by us using Schaller et al. (1992) and Kurucz (1992) 
models, are compared with observations in figures (a), 
(b), (c), (d) and (e) respectively as a function of T e ff. 
Observational data are from Flower (1977) or John- 
son (1966) as indicated, with filled circles and inverted 
triangles representing data for main-sequence and su- 
pergiants respectively. The computed quantities for 
main-sequence and supergiants are denoted by solid 
and dashed lines respectively. Lyman continuum sur- 
face fluxes are compared with the computations of 
Mas-Hesse & Kunth (1991). 



Fig. 4. — Evolution of (a) H/3 to B band luminos- 
ity ratio (log(^ )), (b) B - V, (c) V - R and (d) 
V — K colors for instantaneous burst and continu- 
ous star formation models are shown. Solid and dot- 
ted lines represent smoothed (over 0.5 Myr) and un- 
smoothed IB models. Dot-dashed and dashed lines 
represent CSF models with constant and exponen- 
tially decreasing (e-folding time scale of 10 Myr) star 
formation rates respectively. Note the initial steady 
phase, a red bump and another steady phase during 
the IB color evolution. The IMF used is the same as 
for Fig. 3. 
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Fig. 5. — log(-^ ) is plotted against B — V, V ~R and 

V — K colors in figs (a), (b) and (c) respectively. The 
dotted line is for the instantaneous burst with thick 
dots at every 0.5 Myr interval. Solid lines identified 
by 0, 3.5, 5.0 and 6.0 are tracks for models with co- 
existing stellar populations from two bursts of equal 
strengths. The labeled numbers denote the age of the 
younger of the two populations in Myr, whose posi- 
tions are marked by the asterix on the IB model. The 
older population age varies along the track between 
and 13 Myr. CSF model with constant star forma- 
tion rate is also plotted (dashed line) for 20 Myr of 
evolution. The IMF used is the same as for Fig. 3. 

Fig. 6.— Evolution of (a) log(^ ), (b) B — V, 
(c) V — R and (d) V — K colors for IB models are 
shown for 3 IMFs. The parameters of IMF used are: 
mi = 1, m u = 120 and slopc= 1 (dotted line), to/ = 1, 
m u = 60 and slope= 2.5 (solid line) and mi = 1, 
m u = 30 and slopc= 3.5 (dashed line). Notice that 
as the massive star fraction increases in the cluster, 
colors become bluer in the early phase and redder 
after 5 Myr. 

Fig. 7. — Same as Fig. 6, but for continuous star 
formation scenario with constant star formation rate. 

Fig. 8.— log(-^- ) is plotted against B—V, V-R and 

V — K colors in Figs (a), (b) and (c) respectively for 
the composite two-burst models for the three IMFs 
used in Figs 6 and 7. Burst parameters are identi- 
cal to Fig. 5. Ages of the younger bursts in Myr are 
denoted close to the tracks. Thick dots are placed ev- 
ery 0.5 Myr starting at 6 Myr. IMF differences stand 
out for composite models when the younger burst is 
younger than 3.5 Myr. Beyond this age differences 
tend to decrease especially for V — K color. 

Fig. 9. — Errors due to statistical fluctuations in- 
trinsic to star formation process on bolometric, H/3, 
B band and K band luminosities are plotted for dif- 
ferent cluster masses. Numbers on the curves rep- 
resent mass of the clusters in units of 6641M , the 
mass containing one star earlier than 04 for the IMF 
with mi = 1, m u = 120 and slope= 2.5. The errors 
are computed as the rms fluctuations on the quanti- 
ties over 100 Monte Carlo simulations. See text for 
details. 

Fig. 10. — Same as Fig. 9, but the errors on log(-^- 
), B — V, V — R and V — K colors are plotted. 



Fig. 11. — Comparison of the cluster + nebular spec- 
trum for an IMF with m; = 2, m u = 120 and 
slope = 2.0 at ages 2, 5, 7.5, 10 and 20 Myr are done 
in Figs (a), (b), (c), (d) and (e) respectively for the 
instantaneous burst model. Thick lines are from our 
computations, with dots corresponding to the wave- 
lengths at which computations are performed. Dotted 
lines are from Mas- Hesse & Kunth (1991). Agreement 
is good in general with differences exceeding 0.1 dex 
only at 7.5 and 10 Myr, at which stage cluster is dom- 
inated by fast evolving red supergiants. 
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